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ABSTRACT 

Milli-arcsecond resolution Very Long Baseline Array (VLBA) observations of the archetype WR+0 
star colliding-wind binary (CWB) system WR 140 are presented for 23 epochs between orbital phases 
0.74 and 0.97. At 8.4 GHz, the emission in the wind-collision region (WCR) is clearly resolved as a 
bow-shaped arc that rotates as the orbit progresses. We interpret this rotation as due to the O star 
moving from SE to approximately E of the WR star, which leads to solutions for the orbit inclination 
of 122° ± 5°, the longitude of the ascending node of 353° ± 3°, and an orbit semi-major axis of 
9.0 ± 0.5 mas. The distance to WR140 is determined to be 1.85 ± 0.16 kpc, which requires the O 
star to be a supergiant. The inclination implies the mass of the WR and O star to be 20 ± 4 and 
54±10 Mq respectively. We determine a wind-momentum ratio of 0.22, with an expected half-opening 
angle for the WCR of 63°, consistent with 65° ± 10° derived from the VLBA observations. Total flux 
measurements from Very Large Array (VLA) observations show the radio emission from WR 140 is 
very closely the same from one orbit to the next, pointing strongly toward emission, absorption and 
cooling mechanism(s) that are controlled largely by the orbital motion. The synchrotron spectra 
evolve dramatically through the orbital phases observed, exhibiting both optically thin and optically 
thick emission. We discuss a number of absorption and cooling mechanisms that may determine the 
evolution of the synchrotron spectrum with orbital phase. 

Subject headings: stars:binaries:general — stars:early-type — stars:individual(WR 140) — stars: Wolf- 
Rayet — radio continuum:stars — astrometry 



1. INTRODUCTION 

Wolf-Rayet (WR) stars have dense stellar winds 
from which free-free continuum emission is observed 
from near-IR to radio wavelengths. A significant 
proportion (~ a quarter) of th ese objects also ex- 
hibit non-thermal radio emission llLeitherer et alJ 119951 
11^97; Chaoman et al. 1999: Dou ghertv fc Winiams y2000': 
[CaDDa ct al. 2004), which requires the presence of a pop- 
ulation of rclativistic electrons, in addition to a mag- 
netic field. The acceleration of particles in these ob- 
jects is generally attributed to first-order Fermi ac- 
celeration in shocks within the stellar winds, aris- 
ing either from wind instabilities (..Chen fc White 19941 
or, in the case of binary syst ems, collisions of th e 
stellar winds of the two stars ijEichler fc Usovl Il993() . 
The colliding-wind binary (CWB) model is strongly 

Electronic address: |sean. dougherty @nrc . caj 
^ National Research Council, D.R.A.O 



supported by high-resolution radio, optical and in- 
frared imaging of the very wide WR-I-OB binaries 
WR 146 and WR 147, where the non-thermal emission 
is shown to arise between the WR and OB stars at 
the point of ram-pressure balance between th e two stel- 
lar winds llDqughertv et alJ 119961: IW iUiams e t"anil997t 
iNiemela et al.lll998HDoughertv et alJl20001 

The archetype of CWB systems is the 7.9-year pe- 
riod WR-HO binary system WR 140 (HD 193793). This 
binary comprises a WC7 star and an 04-5 star in a 
highly elliptical orbit (e w 0.88), where the stellar sepa- 
ration varies between ~ 2 AU at periastron to ~ 30 AU 
at apastron. This highly eccentric orbit clearly mod- 
ulates the dramatic variations in the emission from 
the system, observed from X-ray to radio wavelengths 
(see iWilliams et al.l[T990l) . The X -ray luminosity rises 
steadily toward periastron passage, while spectra taken 
just after periastron clearly sh ow extra absorption that is 
visible until about phase 0.1 l|Pollock et al.ll2004^ . The 
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near-IR emission is characterized by a sudden peak in 
emission shortly after periastron, foUowed by a slower 
decline. This behavior is attributed to dust formation in 
the WCR at some point shor tly aft er periastron, and its 
subsequent cooling ( WiUiam^ l2002(l . 

Perhaps the most dramatic variations are observed at 
radio wavelengths (Fig. QJ, where there is a slow rise 
from a low, apparently thermal state close to periastron 
of a few mjy, to a frequency-dependent peak in emis- 
sion of several 10s of mJy between orbital phase (0) 0.65 
to 0.85, before a precipitous decline before periastron 
llWilliams et aLllToOfilWhite fc Beckei<IT995l) . The radio 
variations have been widely attributed to an underlying 
synchrotron source viewed through the changing free-free 
opacity of the extended stellar wind envelopes of the bi- 
nary system alon g the line-of-sight to t he WCR as the 
orbit progresses llWilliams et all 119901: lEichler & Usovl 
IT9931 IWhite fc Beckeilll995D . However, none of the free- 
free opacity models explain the radio light curve in a sat- 
isfactory manner. Alternatively, the variations may be in 
part due to processes intrinsic to the WCR. Changes in 
the stellar separation (D) alter the intrinsic synchrotron 
luminosity of the WCR (oc _D~^/^), and the free-free 
absorp tion and synchrotron se lf-absorption within the 
WCR (|Doughertv et al.ll2003a|) . Furthermore, the im- 
pact of the Razin effect and Coulombic cooling increases 
as the separation decreases. These processes are domi- 
nant at lower frequencies. Inverse-Compton (IC) cooling 
of the shocked gas in the WCR due to the intense ultra- 
violet radiation from the stellar components is important 
at higher frequencies, and also varies strongly with sep- 
aration. 

We report on high resolution radio observations of 
WR 140 obtained with the VLBA that image the WCR 
at a linear resolution of a few AU, approximately the 
stellar separation at periastron. These observations are 
used to constrain the orbit of WR 140, allowing the first 
reliable derivation of the orbital inclination, the longi- 
tude of periastron, and a robust distance estimate. In 
addition to supporting radiometry from the VLA and 
MERLIN, the high resolution observations are essential 
to constrain models of the WCR in WR 140. 

2. OBSERVATIONS 
2.1. VLBA observations 

A 23-epoch campaign of observations of WR 140 at 
both 8421 MHz (X band) and 1667 MHz (L band) was 
carried out using the VLBA (Napier 1995). The cam- 
paign was started on 1999 January 04 near the peak of 
radio emission around orbital phase 0.75. Observations 
continued until 2000 November 18 = 0.97) when the 
radio emission had declined to close to radio minimum 
(see Fig^and could no longer be detected by the VLBA. 

Data were received in dual circular polarizations, and 
acquired from four contiguous 8-MHz bands using four- 
level (2 bit) sampling for each polarization, providing a 
32 MHz total bandwidth. The observations wer e per- 
forme d using phase referencing (see Bcaslev & C onwavl 
1199 5^ where regular observations of a nearby extra- 
galactic background source were made to derive resid- 
ual antenna-based instrumental and atmospheric phase 
errors, which were then interpolated to correct the an- 
tenna gains during observations of WR 140. This tech- 
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Fig. 1. — Radio emission from WR140 at 15 (circles), 5 (trian- 
gles), and 1.5 GHz (squares) as measured with the VLA (solid) and 
WSRT (open) d uring the orbital cycle betwe en 1 984.3 and 1992.2. 
Data t aken from lWilliams et alj tl990l 11994) and IWhite fc Beckeil 
1199^ . 



nique allows the detection of weak emission by enabling 
direct Fourier inversion of the visibility data for imaging 
of WR 140, along with precise astrometric referencing of 
the images. The phase calibrator was J2002-f456, sepa- 
rated from WR 140 by 2.44°. A switching cycle of 90 sec- 
onds on a calibrator and 120 seconds on the target source 
was used, giving ~ 3 hrs on-source time for WR140 at 
each epoch at each frequency. The absolute flux scale 
of the observations was based on internal noise source 
calibration. The data were correlated with the VLBA 
correlator in Socorro, New Mexico, and data calibration 
and imaging was done using standard techniques in the 
NRAO AiPS software package. 

The resulting phase-referenced images at 8.4-GHz for 
each observation epochs are shown in Fig. [21 where the 
rms uncertainty of the image background is between 0.06 
and 0.22 mJy beam~^. The nominal resolution (FWHM) 
of the synthesized beam is 2.0 x 1.5 mas^, which at a 
distance of 1.85 kpc (see §3.2) gives a linear resolution of 
3-4 AU, much smaller than the separation of the stellar 
components in WR 140 for the bulk of the orbit except 
close to periastron. 

The 8.4-GHz emission detected by the VLBA is clearly 
resolved. We identify this emission as arising from the 
WCR in WR 140 since this is the only emission in the sys- 
tem with sufficient brightness temperature (^ 10^"^ K) 
to be detected by the VLBA. The stellar winds, with 
brightness of ^ 10"* K are undetected. A bow-shaped 
ridge of emission is observed at most epochs, as antic- 
ipated for the WCR from model calculations (e.g. see 
'Eichler fc Usov 1993 HCanto et al .'1996: Doug herty et alJ 
■2003a) . with the bow shock wrapping around the star 
with the lower wind momentum - typically the O star. A 
bow-shaped WCR has been suggested from observations 
of other C WB systems (e.g. see [Williams ct al. 199'^ 
iDoughe rtv 2002), but not as clearly as in the observa- 
tions presented here. 

Between orbit phase 0.74 and 0.95, the WCR exhibits 
proper motion that is largely from east to west, moving 
~ 10 mas over the observing period, and also appears to 
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Fig. 2.— 8.4-GHz VLBA observations of WR140 from orbital phase 0.737 through 0.889. The contour levels are -1, 1, 1.6, 2.6, 4.1, 
6.5, 10.5p where p = 220 //Jy beam~^. Note that p is a common contour scale factor between all the epochs to enable easy comparison of 
each image, and is not the rms background level which is between 0.06 and 0.22 mjy beam~^. The greyscale extends from -0.4 (white) to 
2 mJy beam~^ (black). The contour levels and greyscale are identical in each image. The typical beam size in the images is 2.0 x 1.5 mas^. 
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Fig. 2. — (continued) 8.4-GHz VLBA observations of WR 140 from orbital phase 0.900 through 0.973. Contour levels and grey scales are 
the same as Fig. 1^. 
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Fig. 3.— Examples of 1.7-GHz VLBA observations of WR140 at orbit phases 0.806, 0.865, and 0.920. The contour levels are -1, 1, 1.6, 
2.6, 4.1, 6.5p where p = 0.9 mjy. The greyscale extends from -0.1 (white) to 8.3 mjy beam~^ (black). The contour levels and greyscale 
are identical in each image. The typical beam size in the images is 10 X 5 mas^. 
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TABLE 1 

Positions and fluxes from the VLBA observations 



Obs. Code 


Obs. Date 


Phase^ 




8 


.4 GHz data 




1.7 GHz data 












Flux 


PA'^ 








Flux 








(ms) 


(mas) 


Tin .Tv t 


\ ) 


(°) 




(lUcLs) 




1 no A 

DID LKJAA 


1 ooo /m /fi/i 
iyyy/Ui/U4 


V. to t 


U 


U 


10.5 


47 ±4 


( O 


-0.889 


4.45 


12.4 


BB102B 


1999/03/03 


0.757 


-0.003 


0.01 


10.2 


47 ±4 


77 


-0.423 


-1.33 


8.3 


BBllOA 


1999/04/08 


0.769 


-0.045 


-0.34 


10.3 


56 ± 5 


70 


-0.428 


0.86 


14.0 


BB109A 


1999/07/05 


0.800 


-0.106 


0.46 


10.4 


53 ± 4 


71 


-0.526 


0.37 


10.4 


BB109B 


1999/07/23 


0.806 


-0.126 


0.19 


8.3 


48 ± 15 


63 


-0.394 


-0.56 


14.6 


BB109C 


1999/08/12 


0.813 


-0.188 


0.77 


10.2 


53 ± 11 


66 


-0.528 


1.84 


16.8 


BB109D 


1999/09/06 


0.821 


-0.184 


1.44 


11.0 


47 ± 15 


85 


-0.601 


2.70 


18.0 


BB109E 


1999/10/02 


0.830 


-0.232 


0.88 


8.8 


50 ± 15 


77 


-0.699 


0.32 


16.3 


BB109F 


1999/11/01 


0.841 


-0.400 


-0.26 


8.2 


54 ± 10 


73 


-0.773 


0.23 


9.5 


BB117A 


1999/11/21 


0.848 


-0.416 


-0.14 


8.8 


55 ± 11 


88 








BB117B 


1999/12/22 


0.858 


-0.464 


-0.30 


10.5 


55 ±4 


68 


-1.215 


4.75 


16.5 


BB117C 


2000/01/11 


0.865 


-0.450 


-0.42 


9.3 


59 ± 4 


74 


-1.374 


2.96 


17.3 


BB117D 


2000/02/05 


0.874 


-0.469 


-0.73 


8.1 


61 ± 5 


70 






f 


BB117E 


2000/02/27 


0.881 


-0.505 


-0.76 


9.3 


71 ±6 


70 






f 


BB117F 


2000/03/19 


0.889 


-0.581 


0.02 


7.4 


103 ± 15 


90 






f 


BB124A 


2000/04/20 


0.900 


-0.569 


-0.63 


6.6 


70 ±8 


76 






f 


BB124B 


2000/05/26 


0.912 


-0.591 


-1.01 


5.9 


78 ±5 


81 


-1.038 


-0.86 


8.7 


BB124C 


2000/06/19 


0.920 


-0.543 


0.14 


3.7 


90 ±4 


81 


-0.966 


1.93 


5.8 


BB124D 


2000/07/21 


0.931 


-0.629 


0.15 


3.9 


80 ±7 


79 


-1.012 


2.96 


7.2 


BB125A 


2000/08/20 


0.942 


-0.668 


0.41 


3.3 


90 ±4 


76 


-0.963 


2.07 


3.0 


BB125B 


2000/09/25 


0.954 


-0.804 


-0.56 


1.9 


96 ±5 


83 






< 0.3s 


BB125C 


2000/10/22 


0.963 






< 0.3 










< 0.3 


BB125D 


2000/11/18 


0.973 






< 0.2 










< 0.3 



^Based on orbit period and epoch of periastron from lMarchenlco et all i2nn.'yi 
•^Offset from a=20'' 20™ 27.977444'' S= 43°5l'l6. 29174" 
'^Derived using MPS routine .IMFIT 

^Position angle of the orientation of the emission ar c - se c 42.11 
'^O is the observed opening angle of the WCR - see ^3.21 
^Unable to adequately calibrate antenna phase - see ^2.11 
^3cr upper limit 



rotate from "pointing" NW to W. These observations are 
key to unraveling the orbital motion of WR 140 which we 
discuss more fully in ^'6.1\ 

It appears the WCR has a similar extent (and flux) 
from phase 0.737 to 0.841, and then decreases in size, no- 
tably after phase 0.874, and again after phase 0.931. We 
have estimated the brightness temperature of the WCR 
through modelling of the visibilities at 8.4 GHz using the 
DIFMAP analysis package at orbital phases 0.737, 0.858 
and 0.931. At phase 0.737, the WCR extends out to 
~ 3.5 mas from its centre and has a brightness tem- 
perature of ~ 3 X lO"^ K. This brightness temperature 
is essentially constant at phases 0.858 and 0.931, even 
though the WCR is considerably smaller by phase 0.931, 
only extending ~ 1.5 mas from the peak of emission. 

Between phase 0.80 and 0.84, the arc of emission ap- 
pears to develop two peaks along its ridge, which then 
fade away and the arc is reformed by phase 0.86. Dur- 
ing these phases, it appears that the intensity from the 
stagnation point is reduced. The concurrent VLA ob- 
servations (see ti2.2|l show a marked variation in flux at 
1.5 GHz at this stage of the orbit. This emission is still 
brightening from around 18 mJy at phase 0.80 to a peak 
of 26 mJy at 0.83, before fading to 20 mJy by phase 
0.85. This is followed by the suggestion of another slight 
brightening before the onset of the decline to minimum. 
There is little evidence of any similar brightening in the 
emission at the other observed frequencies. 

A similar two-peaked structure was also apparent in 



MERL IN observations of W R 147, a very wide WR+OB 
CWB l|Watson et al.l 12002(1 . It is difficult to attribute 
the reduced level of emission between the peaks seen in 
WR 147 to higher free-free extinction along that line-of- 
sight relative to the peaks because the line-of-sight ex- 
tinction is very low in such a wide binary system. Al- 
ternatively, such a twin-peak structure may be the re- 
sult of time- variable wind density resulting in changes in 
the post-shock density in the WCR. Certainly, between 
phase 0.80 and 0.84 the 15-GHz VLA observations of 
WR 140 (see Fig. ^ have large scatter that may be re- 
lated to such structure changes. We draw attention to 
a similar structure seen in the synthetic 5 GHz image 
of WR 147 produced bv lDoughertv et all l|2003ai Figure 
17). However, in this case, the structure is attributed to 
numerical instabilities in the hydrodynamics code used 
to generate the model (J. Pittard, priv. communication). 

By phase 0.963, the WCR is no longer detected by 
the VLBA. This was expected from the radio light curve 
(see Fig. P) since by this stage of the orbit, the radio 
emission at all frequencies has declined to less than a few 
mJy. Based on our experience from these observations, 
we expect to be able to detect the WCR again somewhere 
around phase 0.3-0.4 during the next orbit. 

The WCR is also resolved in the 1.7-GHz observations 
(see Fig.O. As at 8.4 GHz, proper motion of the source 
is suggested as the orbit progresses. However, the proper 
motion at 1.7 GHz does not appear as continuous as at 
8.4 GHz. Our difficulties with the 1.7-GHz data can be 
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TABLE 2 
VLA Observations of WR 140 



Obs. Date Phase'' 


1.5 GHz 


4.9 GHz 


8.4 GHz 


15 GHz 


22 GHz 


Array 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 



1998/06/28 
1999/05/19 
1999/06/13 
1999/07/06 
1999/07/13 
1999/07/31 
1999/10/03 
1999/10/21 
1999/11/06 
1999/11/21 
1999/12/04 
1999/12/27 
1999/12/29 
2000/01/14 
2000/02/23 
2000/03/29 
2000/04/14 
2000/05/11 
2000/06/21 
2000/07/17 
2000/09/27 
2000/10/15 
2000/11/21 



0.671 
0.783 
0.792 
0.800 
0.802 
0.809 
0.831 
0.837 
0.842 
0.848 
0.852 
0.860 
0.861 
0.866 
0.880 
0.892 
0.898 
0.907 
0.921 
0.930 
0.955 
0.961 
0.974 



4.67±0.10'' 
16.94 ± 1.90 
16.99 ±0.06 
18.50 ±0.07 
18.75 ±0.07 
19.42 ±0.07 
26.37 ±0.07 

24.08 ±0.07 
22.23 ±0.09 
21.18 ±0.12 

20.09 ±0.18 
22.54 ±0.20 

21.58 ±0.17 
20.54 ±0.12 
24.00 ±0.26 
22.22 ±0.39 

20.59 ±0.57 
17.29 ±0.28 
13.56 ±0.44 
8.60 ±0.57 
2.35 ±0.74 
1.20 ±0.09 
0.62 ±0.07 



19.24 ±0.09 
25.96 ±0.12 
25.94 ±0.06 
25.94 ±0.06 
26.50 ±0.07 
25.04 ±0.09 
27.40 ± 0.07 
24.52 ±0.07 
24.73 ± 0.08 
23.48 ± 0.09 
21.87 ±0.09 
22.29 ±0.09 
20.92 ±0.09 
19.99 ±0.06 
20.63 ±0.06 

18.50 ±0.06 

17.51 ±0.06 
15.27 ±0.07 
13.59 ±0.07 
10.00 ±0.07 
4.74 ±0.12 
3.51 ± 0.07 
1.66 ± 0.05 



25.14 ±0.08 
22.80 ±0.07 
22.44 ± 0.06 
21.77 ±0.16 
22.60 ±0.08 
21.49 ±0.08 
22.66 ±0.07 

20.79 ± 0.04 
21.00 ±0.05 

20.38 ±0.06 
18.65 ±0.07 
18.68 ±0.12 

17.80 ±0.09 
16.74 ±0.06 
17.55 ±0.07 
15.76 ±0.07 
15.22 ±0.06 
13.52 ±0.06 

12.39 ±0.09 
9.71 ±0.07 
6.24 ±0.06 
4.48 ± 0.07 
2.64 ±0.06 



19.36 
19.51 
18.78 
20.32 
18.77 
18.61 
18.40 
16.85 
16.17 
17.18 
16.77 



±0.24 
±0.18 
±0.17 
±0.22 
±0.18 
±0.19 
±0.19 
±0.20 
±0.16 
±0.21 
±0.29 



17.91 ± 0.39 



14.57 ±0.23 
13.44 ±0.17 



11.71 ±0.20 
11.81 ±0.20 
9.95 ±0.11 
6.94 ±0.09 



17.87 
18.55 
16.77 
15.20 
16.28 
14.84 
14.65 
16.87 



± 0.42 
± 0.58 
± 0.62 
± 0.52 
± 0.22 
± 0.23 
± 0.21 
± 0.26 



13.46 ± 0.21 
12.59 ± 0.18 

11.79 ±0.24 
11.68 ± 0.21 
9.52 ± 0.16 
10.00 ±0.38 
10.19 ± 0.23 
7.48 ± 0.14 
7.16 ±0.25 
5.00 ±0.11 



BnA 

D 
DnA 

A 

A 

A 
BnA 
BnA 

B 

B 

B 

B 

B 

B 
CnB 

C 

C 

C 
DnC 
DnC 

D 

A 

A 



^Based on orbit period and epoch of periastron from lMarchenko et al 



la errors measured in outlying image regions. 



attributed to rapidly varying differential ionospheric de- 
lay and phase offsets between the position of the phase- 
reference source and WR140 during the observations. 
Throughout the observing campaign, solar activity was 
close to the maximum of its 11-year cycle, contributing 
to significant and rapid ionospheric variability. 

A summary of the positions and fluxes of the VLBA 
data is given in Tabled The fluxes and positions were 
derived at both 8.4 and 1.7 GHz using the AiPS routine 
JMFIT. As a check, fluxes were also determined using 
TVSTAT, where values within ±5% of those derived with 
JMFIT were found. This difference is approximately the 
uncertainty that is to be expected in the absolute flux 
scale calibration. However, no flux uncertainty is quoted 
in Table Q] since the VLBA reso lves out ~ 40% of the 
emission from WR140 (see Sec. I2.2() . The uncertainty 
in the positions is most likely dominated by ionospheric 
delay and phase offset error, certainly at 1.7 GHz. Since 
delay oc and offset oc they are also likely to 

dominate at 8.4 GHz over the Gaussian source fitting er- 
rors, which are a fraction of a milli-arcsecond in these 
images. For epochs where the WCR emission appears 
double-peaked, the mid-point of the two peaks is quoted. 
Also included in Table ^ are our estimates of the posi- 
tion angle of the axis of symmetry of the WCR, along 
with their uncertainty. We also include estimates of the 
opening angle of the WCR. Though not entirely robust, 
we believe each of these estimates are good to first order. 

2.2. VLA observations 

Multi-frequency VLA observations of WR 140 were 
made throughout the VLBA observations to provide con- 
tinuum spectral information. Observations were made 
in standard continuum mode, using two 50-MHz chan- 



nels centered on 1465 (20cm), 4885 (6cm), 8435 (3.6cm), 
14965 (2cm), and 22485 MHz (1.3cm). Antenna phases 
were established via phase-referencing using J2007+404, 
and the absolute flux level from observation of the bright 
source 3C286. Data was calibrated and reduced in the 
standard manner using the AiPS package. Since WR 140 
is unresolved in these observations, fluxes were deter- 
mined using JMFIT assuming a point source. The results 
of the VLA observations are given in Table |3 and dis- 
played in Fig. 21 The quoted uncertainties are the back- 
ground rms of the images. We note that the absolute 
flux scale has an uncertainty of ~ 5%. 

The VLBA fluxes at 8.4 GHz are closely 50% of the 
VLA 8.4-GHz fluxes. Assuming the thermal flux at 8.4 
GHz is that at phase 0.974 (2.64 mJy), we estimate the 
VLBA is resolving out 40% of the synchrotron flux. 
This missing flux is most likely emission arising further 
from the shock apex than the emission detected by the 
VLBA. At 1.7 GHz, the VLBA fluxes have large scatter, 
most certainly due to calibration problems. As a result, 
the relationship to the VLA observations is unclear. It 
appears that 40-50% of the emission is undetected at 
phases around 0.8, though there could be evolution of 
this fraction with orbital phase. At phase 0.75 and near 
phase 0.95, there may be little or no flux resolved out, 
suggesting that at these phases the source may be smaller 
than between these phases. 

In Fig. [SI the new data are directly compared with 
those measured by others from the orbit immediately 
prior to the one presented here. The new data shows 
less scatter than the earlier observations, which we at- 
tribute to improved system performance and better an- 
tenna phase calibration. It is clear from these figures that 
the radio emission from WR 140 very closely repeats from 
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Fig. 4.— The VLA observations of WR 140 from Table |^ as a 
function of orbital phase at 1.5 (squares), 4.9 (triangles), 8.4 (open 
diamonds), 15.0 (circles) and 22 GHz (open pentagons). The data 
points are connected to help guide the eye. Gaps are where data 
were not collected at 15 and 22 GHz. 



one orbit to the other. The radio light curves at each of 
the frequencies are very similar, with the peak and the 
decline to minimum also being very similar. There may 
be an offset between the 1.5 GHz data from phase 0.88 to 
0.92. However, we caution against over interpreting this 
apparent difference. The scatter in the previous 1.5 GHz 
fluxes between orbital phase 0.75 and 0.82 leads us to 
think that the apparent offset may not be significant. 
The close similarity between the fluxes during successive 
orbits at all three frequencies leads us to conclude that 
whatever process (es) govern the nature of the radio emis- 
sion from WR 140, they repeat from one orbit to another, 
suggesting that they are controlled primarily by orbital 
motion. 

2.3. MERLIN observations 

In conjunction with the VLBA observations, observa- 
tions using the Multi-Element Radio-Linked Interferom- 
eter Network (MERLIN) in the UK were obtained in De- 
cember 1999 and February 2000 at 1.66 GHz. The ob- 
servations were approximately 10 hours duration, with 
a total bandwidth of 16 MHz in each of two polariza- 
tions. In addition, observations at 5 GHz on June 1992 
and 1.67 GHz on August 1993 were extracted from the 
MERLIN archive. All of these observations were phase 
referenced using J2007-I-404, and reduced in the stan- 
dard manner with the AiPS program. The absolute flux 
scale was determined through observations of 3C286 and 
bootstrapped to the bright point-sources OQ208 and 
0552-^298. WR 140 is not resolved in the MERLIN ob- 
servations at either 1.6 or 5 GHz. The salient results of 
these observations are summarised in Table 13 where the 
quoted fluxes are determined assuming a point source, 
and taking the rms background level of the images to be 
the flux uncertainty. As in the VLA observations, the 
absolute flux scale has an uncertainty of ~ 5%. 

Comparison of the MERLIN fluxes and the VLA ob- 
servations show very similar flux levels, as expected for 
a source that is unresolved by both of these interferome- 
ter arrays. Since WR 140 is unresolved by MERLIN, we 
estimate the extent of the emission region is less than 
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Fig. 5. — The new VLA observations (solid circles) compared 
with observations from previous orbits at 15 (top), 4.9 (middle) 
and 1 .5 GHz (bottom). Open circles are previous VLA observa- 
tions IWhite fc Becker,1995') and open squares are from the WSRT 
I Williams et alj|199(l 119911 



13 mas at 5 GHz. 

3. THE ORBIT AND DISTANCE OF WR140 

3.1. Orbital parameters of WR 1^0 

Many of the orbital parameters in WR 140, in par- 
ticular the orbital period (P), epoch of periastron pas- 
sage (To), eccentricity (e) and the argument of perias- 
tron (w) are well-determined from radial velocity mea- 
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TABLE 3 
MERLIN Observations of WR 140 



Obs. Date 


Phase" 


1.67 GHz 


5 GHz 






(mJy) 


(mJy) 


1992/06/28 


0.906 




13.2 ±0.2 


1993/08/06 


0.055 


< 0.3 




1999/12/29 


0.861 


20.3 ±0.1 




2000/02/02 


0.873 


24.4 ±0.2 




2000/02/05 


0.874 


24.2 ±0.1 





"Ba sed on or 
1200: 



bit parameters from lMarchenko et all 




80 100 120 

Inclination Angle (degrees) 



Fig. 6. — Solutions for the longitude of the ascending node (solid 
line) and orbit semi-major axis (dotted line) as a function of orbit 
inclination, derived from an observed separation of 12.9 mas at 
a position angle of 151.7° on June 17, 2003 (orbit phase 0.297) 
|Monnier et al. 2004). The uncertainty in the IOTA position angle 
gives a spread in the possible values of Q of closely ±1°, whereas 
the uncertainty in the separation gives ±0.3 mas for the semi-major 
axis uncertainty. 



surements (see iMarchenko et all 120031 and references 
therein). However, the orbital incHnation (i), semi-major 
axis (a) and the longitude of the ascending node (fi) re- 
quire that the system is resolved into a "visual" binary. 
Recently, the two stellar components in WR140 have 
been resolved using the Infrared- Optical Telescope Array 
(IOT A) interferometer at a single epoch (Monnicr et al.l 
|2003). This single observation allows for the first time 
the determination of families of solutions for (i,a, 51). 
However, until further IOTA observations are available, 
the VLBA observations of the WCR offer the only means 
to determine the orbit direction and inclination, and 
h ence 57 and o. 

iMonnier et all l|2004D observed WR140 on June 17, 
2003 to have a separation of 12.9]'Iq'4 "las at a posi- 
tion angle of ISl.T^iJ'g degrees east of north. Assum- 
ing P = 2899 da ys, 'r„ = 244 6147.4, e = 0.881 and 
ll> — 47° (Marchen ko et al1l2003(l . this observation at or- 
bital phase 0.297 gives families of solutions for {i, 17, a) as 
shown in Fig. for 0° < i < 180°. We examine the full 
range of possible inclination angles since the direction of 
relative orbital motion dictates whether i is greater or 
smaller than 90°. 

Low orbit inclinations angles i < 30° (or equivalently 
i > 150°) can be easily ruled out since radial veloc- 




0.75 0.8 0.85 0.9 0.95 

Orbital Phase 

Fig. 7. — The change in the position angle of the axis of symmetry 
of the WCR as a function of orbital phase, as estimated from the 
VLBA 8.4-GHz observations (see Table The solid line is the 
weighted best-fit curve of the position angle of the line-of-centres 
of the two stellar components projected on the plane of the sky 
as a function of orbital phase, deduced using the orbit parameters 
from Marchenko et al. 12003) and Mo nnier et al.. L2004) . The hne 
corresponds to i = 122° and Q = 353°. 



ity variations are observed and low inclinations would 
require the stellar components to have untenably large 
masses. Ultra-violet line observations with lUE, thought 
to have been obtained close to co njunction estimate 
i ~ 38° l|Setia Gunawan et al.l[200l. The same anal- 
ysis with orbital elements from Marchenko et al.l l|2003fl 
give higher estimates of i, with larger uncertainty (P.M. 
Williams, priv. communication). ^archcnko ct aQ 
(2003) estimated 5 ± 15° from their opt ical line study, 
and similar work bv lVarricatt et ahl l)2004() using the He i 
emission line in the near-IR gives an estimate of i ~ 65°. 
Unfortunately, these estimates are difficult to constrain 
well, and all have large uncertainty. 

The VLBA observations of the WCR permit us to con- 
strain acceptable solutions for i. Under the assumption 
that the free-free opacity along the line-of-sight to the 
WCR is sufficiently low as to not impact the apparent 
distribution of emission from the WCR, we expect the 
arc of WCR emission to wrap around the star with the 
lower wind momentum (the O star) and "point" toward 
the WR star. In this case, the rotation of the orientation 
of the WCR as the orbit progresses implies that the O 
star moves from SE to close to due E of the WR star over 
the period of the VLBA observations. If it is assumed 
the axis of symmetry of the WCR emission is coincident 
with the projection on the plane of the sky of the line-of- 
centres of the two stars in the binary, then we can derive 
the orbital inclination from the change in the orientation 
of the WCR with orbital phase. Each {i, 17) family pro- 
vides a unique set of position angles for the projected 
line-of-centres as a function of orbital phase. Using a 
weighted minimum measure to fit the position angle 
of the line of symmetry of the WCR as a function of 
orbit phase for different sets of (i,17), we find a best- 
fit solution of i = 122° ± 5 and 17 = 353 ± 3 (Fig. d). 
These values lead to a value for the semi-major axis of 
a — 9.0 ±0.5 mas, where the uncertainty is dominated by 
the uncertainty in the separation as observed by IOTA. 
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TABLE 4 
Basic parameters of WR140 



20''20"'27, 97800' 37,97750" 27 97700' 27.97650' 

RighL Ascension (J3000) 

Fig. 8. — The orbit of WR140 on the plane of the sky using 
e = 0.88, = 47°, O = 353°, i = 122° and a = 9.0 mas at 
orbital phase 0.737, 0.858 and 0.931, overlaid on the VLBA 8.4- 
GHz images. The WR star is to the W (right) of the WCR at these 
phases. The proper motion and the rotation of the WCR as the 
orbit progresses is clear. The relative position of the stars and the 
WCR was determined using a wind momentum ratio of 0.22 (see 

El- 



Using the semi-major amplitudes of the radial veloc- 
ity cu rve, eccentricity and period from iMarchenko et al.l 
l)2003j) . an orbit inclination of z = 122° ± 5° gives masses 
of 20 ± 4 Mq and 54 ± 10 Mq for the WC and O star re- 
spectively, within the anticipated mass ranges of a WC7 
and an 04-5 star. We note that the uncertainty in the 
masses is dominated by that in the orbital inclination. 

The derived orbit inclination is consistent with values 
previously suggested. However, this newly derived orbit 



Parameter Primary Secondary System 



Distance (kpc) 
M„ 

Spectral type 

BC 

Mboi 

log (Lboi/L©) 
Mass (M0) 

Vcx> (km s-1) 

M (Me yr-1 



-6.4 
04-5 I 
-4.3 
-10.7 
6.18 ° 
54± 10 
3100 
8.7 X 10" 
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-5.6 
WC7 
-3.4 ^ 
-9.0 
5.5 ° 
20 ±4 
2860 ° 
4.3 X 10-5 



1.85 
-6.8 



^Based on and the calibration 
Williams ct al.> (I1990D 
^Calculated assuming MQ^bol — 4.75 
iSetia Gunawan j5^_aLj 



lEencn s & Will iamsH1994ll 



of lVacca et al.l 1199611 
lAitoll976ll 



presents a challenge to current models of dust forma- 
tion. To date, most models of dust formation in WR 140 
assume that the gas from which dust is formed in the 
WCR, is compressed within ~ 0.15 yr of periastron pas- 
sage, when densities in the WCR are sufEciently high 
for dust formation to occur ("WiUiam s et alJ ll990l. The 
subsequent motion of the compressed gas is determined 
by the velocity of this material when it is compressed. 
Since the momentum of the WR star wind is higher than 
that of the O star, this material moves away from the 
WR star, along the WCR. With the orbit orientation 
derived here, the O star is NW of the WR star during 
periastron, and material compressed at periastron will 
therefore have a proper motion to the NW. However, 
this is not what is observed. High-resolution IR observa- 
tions show that dust ejected during the 2001 periastron 
passage has prope r motion from we st to east, away from 
the WR star (Mo nnier et al.ll2002D . New dust models 
are now attempting to address this challenge (Williams 
2005, in preparation). 

If the WR star was to NW of the O star at periastron, 
as previously hypothesized from the dust proper motion, 
the WCR would then have to be wrapped around the 
WR star in the VLBA observations. This requires the 
unlikely scenario of the 0-star wind momentum being 
greater than that of the WR star wind. We acknowl- 
edge that free-free opacity of the stellar winds can effect 
the observed shape of the WCR, but we anticipate the 
free-free opacity in WR 140 would only be sufficiently 
large to have a major impact on the apparent distri- 
bution of the 8.4-GHz emission close to periastron (see 
Dougherty ct al. 2003a, Fig. 11 for an example). How- 
ever, until full hydrodynamical simulations of WR 140 
are available, the full impact of free- free opacity on the 
appearance of the WCR remains unclear. High resolu- 
tion near-IR observations at several bands should be able 
to distinguish the identity of the two stars since we an- 
ticipate the WR star is significantly brighter at longer IR 
wavelengths than the O star. 

3.2. Distance of WRI4O 

Distance estimates of WR stars are typically based 
on absolut e magnitude calibrati ons that often have large 
scatter (see lvan der Huchtl200lD . Having determined the 
orbital inclination and semi-major axis it is now possible 
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to make an estimate of the distance to WR140 indepen - 
dent of any stellar parameters. iMarchcn ko et al.l l|2003f) 
determined a s'mi = 14.10 ± 0.54 AU from radial veloc- 
ity observations, which leads to a = 16.6 ± 1.1 AU for 
i — 122° ± 5°. We have derived the semi-major axis to 
be a = 9.0 ± 0.5 mas, and together these give a distance 
of 1.85 ±0.16 kpc. 

3.3. Basic system parameters of WR 1^0 

This distance is somewhat lar ger than the usually 
quoted value of 1.3 kpc deduced bv lWilliams et alJ ljl99(]D 
from the luminosity of the system. Since the primary 
0-star lumino sity indi cator is masked by the WC7 spec- 
trum, Williams et al.l l|l99Q) assumed a main sequence 
04-5 star with an absolute magnitude of —5.6 and took 
that of the WC7 star to be —4.8. However, they re- 
marked the distance estimate would increase if the O 
star was more luminous than main sequence. With the 
system at 1.85 kpc, the absolute magnitudes of the 04-5 
star and the WC7 star become —6.4 and —5.6 respec- 
tively. An absolute magnitude of —5.6 for the WC7 star 
is bright compared to the mean values (between —4.5 and 
—4. 9) dedu ced from previous estimates for WC7 stars 
(see Ivan dev Hucht 2001, Table 27) but the dispersion is 
large (~ 1™), and certainly the value presented here is 
not extreme. The brightn ess of the 04-5 st ar suggests 
that it is a supergiant fsee lVacca et al.lll996l Table 7). 

With the increase in distance, a reassessment of the 
mass-loss rates of the two stars is appropriate. To cir- 
cumvent any uncertainties presented by clumping cor- 
rections that are required when using radio flux de- 
terminations of mass loss, we turn to X-ray observa- 
tions in order to try and determine the mass-loss rate 
of the WR star. The X-ray luminosity in WR 140 is 
domin ated by emission from th e shocked WR wind ma- 
terial llPittard StevenFil2002() . where any clumps are 
assumed to have b een d estroyed by the shock. From 
IZhekov fc Sk"hl^ (|2003), the X -ray luminosity mea- 
sured by ASCA gives a mass-loss rate for the WR 
star at 1.85 kpc of 4.3 x 10"^ M© yr"!. This can 
be readily reconciled with the un clumped estimate of 
5.3 X 10-5 Mq yr^i from Willia ms etaP ((1990) with 
a reasonable filling factor. Since mass-loss rates derived 
from thermal free- free e missio n cx D^'^, the unclumped 
value of fWiUiams et all l|199(]n for WR 140 at 1.85 kpc 
becomes 9.7 x 10-^ M© yr"^. Comparison with the 
clumping-free. X-ray der ived value implies a fi lling fac- 
tor ~ 0.2. For the O star. lRepolust et all l)2004h suggests 
values of 8.6 — 8.8 x 10"^ M© yr~^ for 04-5 supergiants. 

These mass-loss rates imply a wind momentum ratio 
rj = 0.22, allowing us to position the stars relative to 
the WCR in the VLB A images (see Fig.lHJ. This wind- 
momentum ratio is considerably higher than the 0.035 
deduced by Williams et al. (1990), assuming a main se- 
quence O star, and a mass-loss rate for t he WR star that 
was uncorrected for clumping. Following lEichler fc Usovl 
(21)93), 77 — 0.22 implies a half-opening angle for the 
WCR of 63°. Using the VLB A images in Fig. H we can 
estimate the half-opening angle 6 from 




sin'(/'tan(e/2). 



where ip is the angle between the line-of-sight and the 
line-of-centres between the two stars, and O is the ob- 
served full opening angle of the WCR. Using the orbit 



Fig. 9. — The angle between the Une-of-sight and the Une-of- 
centres of the two stars in WR 140 as a function of orbital phase. 
The vertical lines at phases 0.24 and 0.99 are when this angle tran- 
sits through 63° , the half-opening an gle of the WC R shock cone as 
deduced from r/ = 0.22 and following lEichler fc Usov. l,1993i) . 



derived in ii3.1l to derive ip at the observed orbital phases 
(Fig. EJ, and the values of O given in Table we find 
9 — 65° ± 10°, consistent with the estimate derived from 
the wind- momentum ra tio. Adding further support to 
this value, iVarricatt et a l. (2004) would derive a similar 
opening angle from interpretation of He I line observa- 
tions using the orbital inclination derived here. 

With the orbit size and wind-momentum ratio we 
can compare the observed size of the emission region 
wit h Trrp, the size of th e emission region derived in 
the lEichler fc U"swl l|1993D CWB model. For a stellar 
separation of 26.45 AU at phase 0.737, and 77 = 0.22, 
nro = 26.6AU. The size of the emission region (l) 
and its projection on the sky is simply related by 
/ = l±sec'4!, where ip is defined above. At phase 0.737, 
the size observed by the VLB A is 6 mas, which at 1.85 
kpc corresponds to 11.1 AU. Since ip = 43° at this phase, 
/ ~ 15.2 AU, considerably smaller than expected. To at- 
tain the predicted size would require the observed size 
to be 10.5 mas. As discussed earlier, the VLBA does re- 
solve out ~ 40% of the flux from the WCR in WR 140, 
so the observed size is expected to be smaller than that 
predicted. However, the MERLIN observation imposes 
an upper limit of 13 mas, consistent with the predicted 
size. 

4. ASTROMETRY OF THE WCR 

4.1. Positions and proper motion of the WCR 

The use of phase-referencing permits accurate astrom- 
etry of the radio emission of WR140. The VLBA po- 
sitions for the WCR are given in Table Q The offsets 
of the WCR in the 8.4 GHz observations from the po- 
sition at the first observation epoch as a function of 
time are shown in Fig. 110! A largely westward proper 
motion in the WCR is evident in the 8.4-GHz data. 
A linear regression fit to the 8.4-GHz positions gives 
a proper motion of /i^ = —5.20 ± 0.27 mas yr~^ and 
fig — —0.38 ± 0.27 mas yr~^. In an attempt to assess 
if the ionosphere has any clear impact on the 8.4-GHz 
astrometry, observations between 1800 and 0600 MST 
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Fig. 10. — The positions of tlic WCR at 8.4 GHz as a function 
of time for motion east-to-west (top) and north-to-soutli (bottom). 
Both plots have the same scales for easy comparison. Circled points 
represent observations made between 1800-0600 MST. Solid lines 
are regression fits to all the data points, giving /la = —5.20 it 
0.27 mas yr~^ and /ig = —0.38 ± 0.27 mas yr~^. 



(Socorro) are identified. The east-west positions derived 
from night observations appear to be eastward of the 
daytime positions. Similarly, the north-south positions 
measured during the night are typically northward of day 
time positions. We have no clear explanation for this 
apparent offset in day versus night observations. How- 
ever, repeating the regression fits using only observations 
obtained during the night gives a proper motion fia = 
-5.44 ±0.25 mas yr"! and fj.s = -0.84 ±0.46 mas yr^^. 
The east-west proper motion is closely the same as that 
derived from all data points, though the derived north- 
south motion is larger, albeit with a large uncertainty. 

Hipparcos observed a proper motion in WR 140 of 
Ha = —5.36 ± 0.58 mas yr^^ and us = —2.37 ± 
0.49 mas vr~^ b etween JD 2447870 and 2449046 
ijPerrvman fc ESAII19971) . This date range spans orbital 
phase 0.26 to 0.67. Using the proper motion deduced 
for the WCR together with the orbit defined here and 
T] — 0.22, the deduced proper motions of the O and 
WR stars are fia — —5.21 and —5.91 mas yr~^ and 
Hs = —0.51 and —1.54 mas yr~^ respectively. The de- 
duced east-west motion of both the WR and O stars are 
consistent with the Hipparcos proper motion. However, 
it would appear the deduced north-south motion of the 
WR is closer than that of the O star to the Hipparcos 
observation. This is at odds with the expectation that 
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Fig. 11. — The east- west proper motion of the WCR at 
8.4 GHz (solid circles) and 1.7 GHz (solid squares). An offset 
between the 1.7 and 8.4-GHz positions is clear. Circled points 
represent observations made between 1800-0600 MST. Solid lines 
are regression fits to all the data points at each band, assuming 
fia = —5.20 mas yr~^, the proper motion derived at 8.4 GHz. The 
offset in the two lines is 5 ± 1 mas. 



the Hipparcos observations reflect the proper motion of 
the O star, since the O star is expected to be the dom- 
inant emitter at visible wavelengths (closely Johnson B 
and V) where the Hipparcos observations were made. 
However, given the large uncertainty in the Hipparcos 
north-south proper motion, arguably both deduced val- 
ues are consistent with the Hipparcos observation within 
the measurement uncertainties. 

4.2. Offset between 8.4 and 1.7-GHz emission? 

In spite of the large scatter in the 1.7-GHz positions, 
the 1.7-GHz emission appears to be westward of that at 
8.4 GHz. In Fig. [ill the 1-7 and 8.4-GHz proper mo- 
tion in the east-west direction is shown. From the figure, 
the offset in right ascension at the two bands is read- 
ily apparent. We also note that the largest excursions 
of the 1.7-GHz positions occur during daytime obser- 
vations, lending credibility to our suggestion that iono- 
spheric effects dominate the uncertainty in the astrome- 
try. Taking /i^ = —5.20 as derived from all the 8.4-GHz 
positions, an offset of 5 ± 1 mas is derived between the 
1.7 and 8.4-GHz right ascension values. If we only use 
nighttime data at both bands along with fia = —5.44, the 
position offset between the two bands is 4 ± 1 mas. We 
conclude that despite the poor astrometry at 1.7 GHz, a 
right ascension offset of 4-5 mas between the peak of the 
1.7 and 8.4-GHz emission is observed. 

5. THE RADIO SPECTRA OF WR 140 

The new VLA data at five frequencies allow us to ob- 
serve the radio spectrum and its evolution better than 
previously possible, most particularly at the higher fre- 
quencies. The radio spectra at a number of orbital phases 
are shown in Fig. 1121 and exhibit considerable evolution 
as the orbit progresses. At phase 0.67, the spectrum 
shows a turnover at 8.4 GHz from a positive spectral 
index^ at low frequencies to a negative index at higher 

Flux Sjj ai , where a is the spectral index 
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Fig. 12. — The observed radio spectra of WR 140 at a number of 
orbital phases as observed by the VLA. The orbital phase of each 
spectrum is shown on the left. Not all observed phases are shown 
for clarity. 



frequencies. This is consistent with the data shown in 
Fig.n where a turnover between 15 and 5 GHz is appar- 
ent around phase 0.65. Prior to this phase, the turnover 
frequency is above 15 GHz. As the orbital phase pro- 
gresses after phase 0.67, the turnover frequency moves to 
lower frequency, eventually to below 1.5 GHz by phase 
0.83, when the spectrum has a negative spectral index 
across the observed frequencies. The flux decreases at 
all frequencies beyond this phase and the negative index 
remains until phase 0.907, where the low frequency spec- 
trum evolves rapidly. By phase 0.930, the spectrum is 
flat, and by 0.955 the turnover is greater than 22 GHz, 
and the spectrum exhibits a positive index. 

At phase 0.974, the spectrum is a power-law with a 
best-fit spectral index of 0.72±0.03, a value characteristic 
of the stellar winds in WR+OB binary systems ( William^ 
Il996fl . Comparison with the observations from the pre- 
vious orbit show that the radio flux at 5 GHz is at its 
minimum value at this phase. The brightness tempera- 
ture at this phase is less than 10^ K (undetected by the 
VLB A), implying thermal emission and we will assume 
the emission at this phase represents the thermal emis- 
sion from the stellar winds alone. Certainly, the fluxes 
at phase 0.974 are consistent with the thermal emission 
expected from the stellar winds of a WR-l-0 binary with 
T~ 10"* K, mass-loss rates and terminal velocities as 
given in Tabled) with a filling factor ~ 0.15 — 0.3, similar 
to that derived earlier. Assuming the thermal emission 
from WR 140 is essentially constant throughout the or- 
bit, the synchrotron spectra at each observed phase can 
be determined by simply subtracting the thermal flux at 
phase 0.974 from the total flux. No 15-GHz observation 
was made at phase 0.9 74 so a value of 3 . 8 mJ y is used, 
derived from the data of lWhite fc Beckeil l)1995(l and con- 
sistent with the phase 0.974 spectrum shown in Fig. 1121 
The resulting synchrotron spectra are shown in Fig. 1131 

The synchrotron spectra between phases 0.67 and 0.92 
is optically thin at several frequencies, with a spectral 
index that appears to be closely constant with a — 
—0.5 ± 0.1, very similar to that expected for first-order 
Fermi acceleration of electrons in strong, non-relativistic 
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Fig. 13. — Spectra of the synchrotron emission in WR 140 at 
several orbital phases, determined by subtracting the thermal spec- 
trum at phase 0.974 (presumed to be solely from the stellar winds) 
from the spectra at each phase. The phase of each spectrum is 
shown on the left. The thick dashed line at the top indicates a 
slope of —0.5. 



shocks (e.g. see IDrurvlll983t IJones fc Ellis on '19911 and 
references therein). The turnover frequency between the 
optically thin and thick parts of the spectra appears to 
be related to orbital phase. Between phase 0.67 and 
0.85 the turnover is between 1.5 and 22 GHz, moving 
to lower frequencies as the orbit progresses. Prior to 
0.64, the turnover frequency is greater than 15 GHz (see 
Fig.nj and after phase 0.85 it is less than 1.4 GHz. After 
phase 0.91, the spectra start to flatten as the luminosity 
decreases. By phase 0.96 the turnover frequency is at 
8.4 GHz. 

The optically thick part of the synchrotron spectrum 
has been widely attributed to free-free absorption in 
the stellar winds along the line-of-sight to the WCR. 
[Williams ct al. (1990) found that if the stellar wind en- 
velope is assumed to be radially symmetric, the free-free 
opacity to the apex of the WCR (i.e. near the stagna- 
tion point of the two colliding winds) only exhibits large 
values and variations close to periastron, when lines-of- 
sight to the apex of the WCR have to traverse the in- 
ner regions of the WR star wind. Their modelling of 
two-frequency data led them to suggest a more complex 
geometry with the 0-star envelope providing a lower den- 
sity region extending away from the WR star, very sim- 
ilar to the now generally accepted geometry for the stel- 
lar wind en velopes in CW Bs (c.f . 'Eichlcr fc Usov Ml993t 
iDoughertv et al. 2003al. IWhiTe fc Becker (1995,) pro- 
posed an ad-hoc model to explain the radio peak near 
phase 0.7 using free-free opacity variations alone, where 
the WR wind is equatorially enhanced and masks the 
0-star and the WCR for the bulk of the orbit. In this 
model, the rise to the radio peak at phase 0.7 indicates 
lower free-free opacity as the 0-star passes through the 
disk toward the observer. 

Unfortunately, these models are too simple to explain 
the radio observations of WR 140, as readily acknowl- 
edged by their authors. The VLBA observations show 
the WCR as a distributed emission region and the lines- 
of-sight to the WCR traverse different regions of the stel- 
lar winds of the two stars. As a result, the emerging 
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emission will be a combination of both optically thick 
and thin emission, since even though lines-of-sight to 
the apex may be optically thick, a substantial amount 
of emission arises from optically thin lines-of-sight to the 
downstream flow. Using the newly derived orbit and as- 
suming a half-opening angle for the WCR of 63°, we now 
know the lines-of-sight to the WCR traverse the 0-star 
wind between orbit phases 0.24 and 0.99 (Fig.|2|) during 
which the most dramatic changes in the radio emission 
are observed. Clearly, if stellar wind free-free opacity is 
important it is the 0-star wind opacity that is of most 
concern, not that of the WR star. 

Another shortcoming of previous models is the assump- 
tion the synchrotron emission is optically thin at all ob- 
serving frequencies, and at all orbital phases. Though the 
high frequency data presented here suggests an optically 
thin component at some orbital phases and frequencies, 
it is necessary to account for the optically thick com- 
ponent of the spectra observed during the bulk of the 
orbit. This may, at least in part, to be due to mech- 
anis ms intrinsic to the W CR. Certainly the models o f 
both lWilhams et alJ Eml and lWhite fc Becked (EM^) 
suggest a synchrotron source that is varying dramatically 
throughout the orbit, variations that are modulated by 
the orbital motion. 

Before exploring the potential impact of a few mecha- 
nisms, it is worthwhile estimating the equipartition mag- 
netic field (Beq) at the apex of the WCR. The syn- 
chrotron luminosity (L) at phase 0.85 based on the spec- 
trum in Fig. El is 1.4 x 10"^° erg s~^, and for a source 
size (l) of 7 mas at 1.85 kpc, B^q ^ 50 mC at the apex 
of the WCR. Since B^q oc L^/^r^/^ (jPacholczvkllTgfnl) . 
our estimate is a lower limit since the luminosity is a 
lower limit (the luminosity will increase when observa- 
tions at frequencies greater than 22 GHz are incorpo- 
rated). By phase 0.93, the source size (3 mas) and lumi- 
nosity (5.7 X 10^^ erg s~-^) give B^q ~ 75 mC. Encour- 
agingly, both these estimates are the same order of mag- 
nitude as those determined by Dougherty ct al. (20031a)- 
It is interesting to note that these magnetic field values 
are closely proportional to 1/D, as expected in a dipole 
field. 

Inverse Compton (IC) cooling has been suggested as 
important in WR 140 , with the close proximity of the 
O star to the WCR l|Doughertv et a l. 2003a). The IC 
cooling time at the apex of the WCR is given by 

t/c = 2.378 X 1022^0 (-)^/^ 
Lo ly 

where tq is the distance from the O star to the WCR, 
Lo the bolometric luminosity of the O star, and B is 
the magnetic field. Using the parameters in Tabled at 
phase 0.85 the IC cooling time is ~ 16 hours at 22 GHz. 
This is short compared to the adiabatic cooling time of 
~ 90 hours. Even at apastron, IC cooling dominates 
{tic ^ 170 hours) over adiabatic cooling 430 hours). 
Since IC cooling oc r'^/^ compared to adiabatic cooling 
(oc r), at some point in the downstream flow adiabatic 
cooling will dominate. Synchrotron cooling is not im- 
portant in WR140, with timescales of > 200 days. The 
signature of IC cooling is a spectrum that is reduced 
at all frequencies, but most severely at high frequencies, 
and is not characterised by a high-frequency knee (Pit- 
tard et al. 2005, in preparation). This process may be 



the underlying cause of the reduction in the luminosity 
as the orbit of WR 140 progresses beyond 0.82, given its 
domination over the other cooling mechanisms. 

The presence of thermal ions in the post-shock flow 
can result in a departure of the refractive index of the 
plasma from unity, giving rise to the Razin effect (see e.g. 
IHornbv fc Williamslll966L and references therein). The 
frequency at which such a departure from unity is sig- 
nificant is given by vji w 20(ne/ B), where ng and B are 
the number density of electrons and the magnetic field 
strength respectively. At phase 0.85 and for r] = 0.22, 
the distance from the WR star to the WCR is 14.5 AU. 
A mass-loss rate of 4.2 x 1O~^M0 yr~^ implies a post- 
shock density in the WCR of 6.2 x 10*^ cm"'''. Taking B 
to be 50 mG, ^i? ~ 2.4 GHz, within our observing range. 

For a field of 50 mG, a source of radius 3 mas and a 
flux of 15 mJy at 22 GHz in the optically thin part of 
the spectrum, synchrotron self-absorption will produce a 
turnover (^s) near 0.3 GHz, outside the range of our ob- 
servations. Since Vs esc B^^^R~'^^^, which is very weakly 
dependent on the magnetic field, the source would have 
to be 20 — 30 X smaller for self-absorption to cause a 
turnover within the observed frequency range. 

In highly eccentric systems such as WR140, the rela- 
tive importance of the various emission, absorption and 
plasma processes changes with stellar separation through 
variations in electron and ion density (oc D~^), mag- 
netic field strength (oc D~^), and radiation field density 
(oc D~^). Given these dependences, the asymmetry of 
the radio emission in WR 140 about periastron/apastron 
remains a puzzle - clearly it cannot be explained by any 
one of these processes alone. The evolution of the spec- 
trum is most likely a result of the combined effects and 
relative changes in these mechanisms as the orbit pro- 
gresses. The full impact of these effects on the observed 
emission from WR 140 awaits more sophisticated models 
than used in previous work. 

New radiative transfer models, based on a fully consis- 
tent hydrodynamic treatment of the WCR, have started 
to explore the impact of a number of processes on the 
radio emission from CWBs, including free- free opacity in 
the stellar winds, synchrotron self-absorption, Coulom- 
bic cooling through interactions with post-shock ions, 
plasma effects such as the Razin effect, and Inverse 
Compton cooling by the intense ultra-violet radiation 
field of the nearby massive stars. These models have been 
very successful in explaining th e radio emission from ver y 
wide CWBs such as WR147 l|Doughertv et alJl2003aj) . 
and are now maturing to the point where they will pro- 
vide more insight to the mechanisms acting in systems 
hke WR140 (Pittard et al. 2005, in preparation). The 
observations presented here represent the key constraints 
for these new models. 
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